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In a temporal network, the presence and activity of nodes and links can change through time. To describe 
temporal networks we introduce the notion of temporal quantities. We define the addition and multiplication 
of temporal quantities in a way that can be used for the definition of addition and multiplication of temporal 
networks. The corresponding algebraic structures are semirings. The usual approach to (data) analysis of tem¬ 
poral networks is to transform it into a sequence of time slices - static networks corresponding to selected time 
intervals and analyze each of them using standard methods to produce a sequence of results. The approach 
proposed in this paper enables us to compute these results directly. We developed fast algorithms for the pro¬ 
posed operations. They are available as an open source Python library TQ (Temporal Quantities) and a program 
lanus. The proposed approach enables us to treat as temporal quantities also other network characteristics such 
as degrees, connectivity components, centrality measures. Pathfinder skeleton, etc. To illustrate the developed 
tools we present some results from the analysis of Franzosi’s violence network and Gorman’s Reuters terror 
news network. 
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I. INTRODUCTION 


In a temporal network, the presence and activity of 
nodes and links can change through time. In the last 
two decades the interest for the analysis of temporal 
networks increased partially motivated by travel-support 
services and the analysis of sequences of interaction 
events (e-mails, news, phone calls, collaboration, etc.). 
The approaches and results were recently surveyed by 
Holme and Saramaki in their paper ini and the book 

d. 

Most of temporal social networks data contain the 
information about activity time intervals of their links, 
sometimes augmented by the activity intensity. The 
usual approach to the (data) analysis of temporal net¬ 
works is to transform it into a sequence of time slices - 
static networks corresponding to selected time intervals 
- see for example ifTfil fT9l |23l . Afterward each time 
slice is analyzed using the standard methods for analy¬ 
sis of static networks. Finally the results are collected 
into a temporal sequence of results. In this paper we 
propose an alternative approach, based on the notion of 
temporal quantity, that bypasses explicit construction of 
time slices. The developed algorithms are transforming 
temporal networks directly into results in the form of 
temporal quantities, vectors, temporal vectors or parti¬ 
tions, and temporal networks. 

In the paper, we first present the basic notions about 
temporal networks. In Section III we introduce the 


temporal quantities and propose an algebraic approach, 
based on semirings, to the analysis of temporal net¬ 


works. In the following sections we show that most 
of the traditional network analysis concepts and algo¬ 
rithms such as degrees, clustering coefficient, closeness, 
betweenness, weak and strong connectivity, PathFinder 
skeleton, etc. can be straightforwardly extended to their 
temporal versions. 


II. DESCRIPTION OF TEMPORAL NETWORKS 


For the description of temporal networks we propose 
an elaborated version of the approach used in Paiek ll24l . 
In our approach we also consider values of links (in most 
cases measuring the intensity/frequency of the activity). 
Pajek supports two types of descriptions of temporal 
networks based on presence and on events (Pajek 0.47, 
July 1999). Here, we will describe only the approach to 
capturing the presence of nodes and links. 

A temporal network Mr = (V, £, T, 7^, W) is ob¬ 
tained by attaching the time, T, to an ordinary network, 
where T is a set of time points, t G T- V is the set of 
nodes, C is the set of links, V is the set of node proper¬ 
ties, and yy is the set of link properties or weights 0. 
The time T is usually either a subset of integers, T C Z, 
or a subset of reals, T C M. In Pajek T C N. In a gen¬ 
eral setting it could be any linearly ordered set. 

In a temporal network, nodes v G V and links I G C 
are not necessarily present or active at all time points. 
Let T{v), T G V,ht the activity set of time points for 
the node u; and T{1), T G W, the activity set of time 
points for the link 1. The following consistency condi¬ 
tion is imposed: If a link l{u,v) is active at the time 
point t then its end-nodes u and v should be active at 
the time t. Formally we express this by 
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T{l{u,v)) C T{u)nT{v). 
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The activity set T(e) of a node/link e is usually 
described as a sequence of activity time intervals 
([■Si, fi))i=i^ where Si is the starting time and is the 
finishing time. 

We denote a network consisting of links and nodes ac¬ 
tive in the time t G Thy Af{t) and call it the (network) 
time slice ox footprint of t. Let T' C T (for example, a 
time interval). The notion of a time slice is extended to 

r'by 

Af(r') = U m- 

teT' 

A. Examples 

Let us look at some examples of temporal networks. 

Citation networks can be obtained from biblio¬ 
graphic data bases such as Web of Science (Knowl¬ 
edge) and Scopus. In a citation network JV = 
{V,C,T,V,yV), its set of nodes V consists of se¬ 
lected works (papers, books, reports, patents, etc.). 
There exists an arc a{u, v) G C iff the work u cites 
the work v. The time set T is usually an interval 
of years in which the works were 

published. The activity set of the work v, T{v), is 
the interval and the activity set 

of the arc a{u,v), T{a), can be set to the interval 
\yearp^^i^{u),yearp^^i^{u)] (instances approach) or to the 
interval [yearyear (cumulative approach). An 
example of a property p G V is the number of pages or 
the number of authors. Other properties, such as work’s 
authors and keywords, are usually represented as two¬ 
mode networks. 

Project collaboration networks are usually based 
on some project data base such as Cordis. The set of 
nodes V consists of participating institutions. There is 
an edge e(u: v) G C iff institutions u and v work on a 
joint project. The time set T is an interval of dates/days 
[dayji^^,,dayi^^f] in which the collaboration data were 
collected. T{v) = T and T(e) = {[■s,/] : there ex¬ 
ists a project P such that u and v are partners on P; s is 
the start and / is the finish date of P}. 

KEDSAVEIS networks are networks registering po¬ 
litical events in critical regions in the world (Middle 
East, Balkans, and West Africa) on the basis of daily 
news. Originally they were collected by KEDS (Kansas 
Event Data System). Currently they are hosted by Pams 
Analytical Systems. The set of nodes V contains the in¬ 
volved actors (states, political groups, international or¬ 
ganizations, etc.). The links are directed and are de¬ 
scribing the events: 

(date, actori, actor 2 , action) 

on a given date the actori made the action on the 
actor 2 - Different actions are determining different re¬ 
lations - we get a multirelational network with a set of 


links partitioned by actions C — {Ca '■ a G Actions}. 
The time set is determined by the observed period T = 
[dayfirst, day last]- Since most of the actors are exist¬ 
ing during all the observed period their node activity 
time sets are T{v) = T. Another option is to con¬ 
sider as their node activity time sets the period of their 
engagement in the region. The activity time set T{1) 
of an arc l{u,v) G Ca contains all dates - intervals 
[day, day + 1) - in which the actor u made an action 
a on the actor v. Another possibility is to base the de¬ 
scription on a single relation network and store the in¬ 
formation about the action a as a structured value in a 
triple {day, day + 1, value) 

value = [{actioni, counti),..., {actiouk, countk)] 

and introduce an appropriate semiring over such values 
(see Sectionjl^. 

There are many other examples of temporal networks 
such as: genealogies, contact networks, networks of 
phone calls, etc. 


III. TEMPORAL QUANTITIES 

Besides the presence/absence of nodes and links also 
their properties can change through time. To describe 
the changes we introduce the notion of a temporal quan¬ 
tity a with the activity set Ta ‘fT 

a'{t) t GTa 
§€ tGT\Ta 

where a'{f) is the value of a at an instant t, and §€ de¬ 
notes the value undefined. 

We assume that the values of temporal quantities be¬ 
long to a set A which is a semiring {A, ©, ©, 0,1) for 
binary operations © : A x A —A and © : A x A —>■ A 
mM- This means that (A, ©, 0) is an Abelian monoid 
- the addition © is associative and commutative, and 
has 0 as its neutral element; and (A, ©, 1) is a monoid - 
the multiplication © is associative and has 1 as its neu¬ 
tral element. Also, multiplication distributes from both 
sides over the addition. Note that 0 and 1 denote the two 
elements of A that satisfy the required properties. In 
expressions the precedence of the multiplication © over 
the addition © is assumed. We can extend both opera¬ 
tions to the set A^ = A U {3€} by requiring that for all 
a G Ajs it holds 

a©§f = 3€©a = a and a©9€ = 3€©a = 3€. 

The structure (Ajj, ©, ©, 1) is also a semiring. 

The “default” semiring is the combinatorial semiring 
(Kjj", + , •, 0,1) where + and • are the usual addition and 
multiplication of real numbers. In some applications 
other semirings are useful. 
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FIG. 1. Semiring addition and multiplication in networks. 


In applications of semirings in the analysis of graphs 
and networks the addition © describes the composition 
of values on parallel walks and the multiplication © de¬ 
scribes the composition of values on sequential walks 
- see Figure For the combinatorial semiring these 
two schemes correspond to basic principles of combina¬ 
torics: the Rule of Sum and the Rule of Product ll2^ . 

The semiring (Mq , min, +, oo, 0), = Mg U {oo}, 

is suitable to deal with the shortest paths problem in net¬ 
works; and the semiring ({0,1}, V, A, 0,1) for reacha¬ 
bility problems. The standard references on semirings 
are Q and iflTll . 


A. Semiring of temporal quantities 

Let denote the set of all temporal quantities 

over in the time T- To extend the operations to net¬ 
works and their matrices we hrst dehne the sum (parallel 
links) a © 6 as 

(a © b)(t) = a(t) © b{t) 

and Ts = U T^; and the product (sequential links) 
a Qb as 

(a © b){t) = a{t) © b{t) 
and Tp = Ta n Tb. 

In these definitions and also in the following text, to 
avoid the ‘pollution’ with many different symbols, we 
use the symbols © and 0 to denote the semiring opera¬ 
tions. The appropriate semiring can be determined from 
the context. For example, in the dehnition of addition of 
temporal quantities the symbol © on the left hand side 
of the equation operates on temporal quantities and the 
symbol © on the right hand side denotes the addition in 
the basic semiring 

Let us define the temporal quantities 0 and 1 with 
requirements 0(f) = i€ and 1(f) = 1 for all f G 
T- It is a routine task to verify that the structure 
{A^{T), ©, ©, 0,1) is also a semiring, and therefore so 
is the set of square matrices of order n over it for the 
addition A © B = S 

^ij — CLij © b^j 


and multiplication A © B = P 

n 

k=l 

Again, the symbols © and © on the left hand side op¬ 
erate on temporal matrices and on the right hand side in 
the semiring of temporal quantities. 

The matrix multiplication is closely related to trav¬ 
eling on networks. Consider an entry pij in an instant 
f 

n 

Pij (^) ^ ® {t) © bkj {t) . 

k^l 

For a value pij (f) to be defined (different from §€) there 
should exist in the instant f at least one node k such that 
both the link (f, k) and the link (fc, j) exist - the transi¬ 
tion from the node i to the node j through a node k is 
possible. Its contribution is aik{t) © bkj{t). This means 
that the matrix multiplication is taking into account only 
the links inside the time slice JV{t). 


B. Operationalization 

In the following we shall limit our discussion to tem¬ 
poral quantities that can be described in the form of 
time-interval/value sequences 

a = ((/*,t^i))?=i 

where f is a time-interval and Vi is a value of a on 
this interval. In general, the intervals can be of dif¬ 
ferent types: 1 - [s^,fi]; 2 - [si,fi); 3 - (si,/i]; 4 - 
{si,fi). Also the value Vi can be structured. For ex¬ 
ample Vi = {wi,Ci,Ti) - weight, capacity and transi¬ 
tion time, or Vi = {di,ni) - the length of geodesics and 
the number of geodesics, etc. We require Si < fi, for 
z = 1,..., A: and Si-i < Si, for z = 2,..., k. 

To simplify the exposition we will assume in the fol¬ 
lowing that all the intervals in our descriptions of tempo¬ 
ral quantities are of type 2 - [si, fi) and /i_i < Si, for 
z = 2,... ,k. Therefore we can describe the temporal 
quantities with sequences of triples 

a = ((si,/i,t^i))J=i- 

In the examples we will also assume that T = 

[iminAmax] © 

To provide a computational support for the proposed 
approach we are developing in Python a library TQ 
(Temporal Quantities). In the examples we will use the 
Python notation for temporal quantities. 

The following are two temporal quantities a and b 
represented in Python as a list of triples 
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ALGORITHM 1. Addition of temporal quantities. 

1: function fo) 

2: if length{a) — 0 then return b 

3: if length{b) — 0 then return a 

4: c^[] 

5: {sa,fa,Va) get{a)\ {sb, fb,Vb) get{b) 

6: while {sa < oo) V (sj, < oo) do 

7: if Sa < Sb then 

8: Sc Sa', Vc <- Va 

9: if Sb < fa then fa t- Sb, Sa -S- Sb 

10: else fa fa', (Sa, fa, Va) get{a) 

11: else if Sa = Sf, then 

12: Sc Sa; fa ^ min(/a, fb) 

13: Va ■(— sAdd{va,Vb) 

14: Sa Sb fa', fd fa 

15: if fd < fb then (sa, fa, Va) ^ get{a) 

16: if fb < fd then {sb, fb,Vb) get{b) 

17: else 

18: Sc Sf,; Va Vb 

19: if Sa < fb then fa -S— Sa', Sb -1— Sa 

20: else/c ^ fb', (sb, fb, Vb) ^ get{b) 

21: c.append{{Sa, fa,Va)) 

22: return standardise) 


a = [ (1, 5, 2) , (6, 8, 1) , (11, 12, 3) , 

(14, 16, 2), (17, 18, 5), (19, 20, 1)] 

b = [ (2, 3, 4) , (4, 7, 3) , (9, 10, 2) , 

(13, 15, 5), (16, 21, 1)] 

The temporal quantity a has on the interval [1, 5) (i.e. 
in instances 1, 2, 3 and 4) value 2; on the interval [6,8) 
value 1; on the interval [11,12) value 3, etc. Outside the 
specified intervals its value is undefined, §€. 

The temporal quantities can also be visualized as it is 
shown for a and b at the top half of Figure]^ 

For the simplified version of temporal quantities we 
wrote procedures sum (Algorithm [^l for the addition 
and prod (Algorithm]^ for the multiplication of tem¬ 
poral quantities over the selected semiring. Because, by 
assumption, the triples in a description of a temporal 
quantity are ordered by their starting times, we can base 
both procedures on the ordered lists merging scheme. 
The basic semiring operations of addition and multipli¬ 
cation are provided by functions sAdd and sMul. 

The function length{a) returns the length (number of 
items) of the list a. The function get{a) returns the cur¬ 
rent item of the list a and moves to the next item; if the 
list is exausted it returns a ‘sentinel’ triple ((X),oo,0). 
The statement [s, f,v) ^ e describes the unpacking of 
the item e into its parts. The statement c.appendie) ap¬ 
pends the item e to the tail of the list c. The function 
standard{a) joins, in the list a, adjacent time intervals 
with the same value into a single interval. 

The following are the sum s and the product p of tem¬ 
poral quantities a and b. They are visually displayed at 
the bottom half of Figure]^ 


ALGORITHM 2. Multiplication of temporal quantities. 

1: function prod(a, 6) 

2: if length{a) ■ length{b) = 0 then return [ ] 

3: C t- [ ]; (sa, fa,Va) ^ gel{a); (sb, fb, Vb) -h- get(b) 

4: while (sa < oo) V {sb < oo) do 

5: if fa < Sb then (sa, /a, Wa) getia) 

6: else if fb < Sa then {sb, fb,Vb) get{b) 

7: else 

8: Sa ^ max(sa, Sb)', fa min(/a, fb)', 

9: Va ■(— sMul(Va,Vb) 

10: c.append{{sa, fa,Va)) 

11: it fa = fa then (Sa, fa, Va) get{a) 

12: it fa = fb then {sb,fb,Vb) get{b) 

13: return standardic) 
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aQb : 
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FIG. 2. Addition and multiplication of temporal quantities. 
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a : 



FIG. 3. Addition and multiplication of temporal quantities - 
growth of size. 


[ (1, 

2, 

2) , 

(2, 3, 

6), 

(3, 

4, 2) 



(4, 

5, 

5) , 

(5, 6, 

3) , 

(6, 

7, 4), 



(7, 

8, 

1) , 

(9, 10, 

2), 

(11 

, 12, 

3) , 


(13, 

14 

, 5) 

, (14, 

15, 

7) , 

(15, 1 

6, 

2) , 

(16, 

17 

, 1) 

, (17, 

18, 

6) , 

(18, 1 

9, 

1) , 

(19, 

20 

, 2) 

, (20, 

21, 

1) ] 




[ (2, 

3, 

8) , 

(4, 5, 

6), 

(6, 

7, 3) 

, 


(14, 

15 

, 10 

), (17, 

18, 

5) , 

(19, 

20, 

1) 


Let la — length{a) and 4 = lengthih). Then, assum¬ 
ing that the semiring operations take constant time each, 
the time complexity of both algorithms is 0{la + h)- 
The example in Figure shows that in extreme cases 
the sum can be almost 4 times longer than each of its ar¬ 
guments, and the product almost twice as long as the 
arguments. If T = [tmin,tmax] C N the length of 
a list describing a temporal quantity can not exceed 


C. The aggregated value 

In some applications over the combinatorial semiring 
we shall use the aggregated value of a temporal quantity 
a = ((si, fi, Vi))^=i- It is defined as 

k 

Sa = ^(/i - Si) ■ Vi 
i=l 

and is computed using the procedure total{a). For ex¬ 
ample So = 23 and E6 = 30. Note that Ea + T,b = 
E{a + b). 

D. Temporal partitions 

The description of temporal partitions has the same 
form as the description of temporal quantities a — 
{{si, fi, Vi))i^i. They differ only in the interpretation 
of values Vi G N. In case of partitions Vi = j means 
that the unit described with a belongs to a class j in the 


time interval [si, fi). We shall use temporal partitions to 
describe connectivity components in Section[X| 

We obtain a more adequate description of temporal 
networks by using vectors of temporal quantities (tem¬ 
poral vectors and temporal partitions) for describing 
properties of nodes and making also link weights into 
temporal quantities. In the current version of the library 
TQ we use a representation of a network Af with its ma¬ 
trix A = [auv] 

_ f w{u,v) {u,v) G C 

otherwise 

where w{u,v) is a temporal weight attached to a link 

{u,v). 

E. Products of a temporal matrix and a temporal vector 

In some applications the product of a temporal matrix 
with a temporal vector is useful. There are two products 

- left and right. 

Let A be a temporal matrix of size n x to, v a vector 
of size n, and u a vector of size to. The product from 
left of A with v, denoted by u = v • A, is defined by 

n 

Uj^^ViQaij, j = 

i=l 

and the product from right of A with u, denoted by v = 
A • u, is defined by 

m 

Vi=^aijQ>Uj, i = 

i=i 

In the TQ library both products are implemented as 
functions MatVecMulL{A, v) and MatVecMulR{A, v). 

If a vector v of size n is considered as a column vector 

- an n X 1 matrix - it holds v • A = (v^ © A)^ and 
A*u = A©u. T denotes the matrix transposition 
operation. 

IV. NODE ACTIVITIES 

In this section we show how we can use the proposed 
operations with temporal quantities (the addition) for a 
simple analysis of temporal networks. 

Assume that the values in temporal quantities Quv 
from a temporal network matrix A are positive real 
numbers measuring the intensity of the activity of the 
node u on the node v. We define the activity of a group 
of nodes Vi on a group V 2 (using the combinatorial 
semiring) as 

act('Fi,V2) = ^ auv 

uGVi VGV 2 
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To illustrate the notion of activity we applied it on Fran- 
zosi’s violence temporal network im. Roberto Fran- 
zosi collected from the journal news in the period Jan¬ 
uary 1919 - December 1922 information about the dif¬ 
ferent types of interactions between political parties and 
other groups of people in Italy. The violence network 
contains only the data about violent actions and counts 
the number of interactions per month. 

We determined the temporal quantities 
pol = act({police}, V) + act(V, {police}), 
fas = act({fascists}, V) + act(V, {fascists}) and 
all = act(V,V). They are presented in Figure 
Comparing the intensity charts of police and fascists 
activity with overall activity we see that most of the 
violent activities in the first two years 1919 and 1920 
were related to the police. In the next two years (1921 
and 1922) they were taken over by the fascists. 


TABLE I. Temporal indegrees and outdegrees for the first ex¬ 
ample network. 



indegrees 


outdegrees 

1 

[ (1, 

9, 

1) ] 

1 

[ (1, 

9, 

1) 

2 

[ (1, 

9, 

2) ] 

2 

[ (1, 

3, 

1) 

3 

[] 




(3, 

9, 

2) 

4 

[ (1, 

3, 

1) , 

3 

[ (1, 

9, 

1) 


(3, 

9, 

2) ] 

4 

[ (1, 

9, 

1) 

5 

[ (1, 

9, 

1) ] 

5 

[ (1, 

5, 

2) 

6 

[ (1, 

9, 

1) ] 


(5, 

9, 

1) 

7 

[ (1, 

5, 

1) , 

6 

[ (1, 

9, 

1) 


(7, 

9, 

1) ] 

7 

[ (1, 

9, 

3) 

8 

[ (1, 

9, 

2) ] 

8 

[ (1, 

9, 

2) 

9 

[ (1, 

9, 

2) ] 

9 

[ (1, 

9, 

2) 

10 

[ (1, 

9, 

3) ] 

10 

[ (1, 

9, 

1) 

11 

[ (1, 

9, 

2) ] 

11 

[ (1, 

7, 

1) 

12 

[] 




(7, 

9, 

2) 

13 

[ (2, 

8, 

2) ] 

12 

[] 



14 

[ (2, 

8, 

2) ] 

13 

[ (2, 

8, 

2) 

15 

[ (2, 

8, 

2) ] 

14 

[ (2, 

8, 

2) 
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[ (2, 

8, 

2) 


V. TEMPORAL DEGREES 

For an ordinary graph with a (binary) adjacency ma¬ 
trix A we can compute the corresponding indegree, i, 
and outdegree, o, vectors using (over the combinatorial 
semiring) the relations 


and 


o = A • e 


where e is a column vector of size n = |V| with all 
its entries equal to 1. The same holds for temporal net¬ 
works. In this case the vector e contains as values the 
temporal unit 1 = [(0, oo, 1)]. 

For a temporal network presented in Figurej^the cor¬ 
responding temporal indegrees and outdegrees are given 
in Table |I] For example, the node 5 has in the time inter¬ 
val [1, 5) outdegree 2. Because the arc (5, 7) disappears 
at the time point 5 the outdegree of the node 5 dimin¬ 
ishes to 1 in the interval [5,9). 

We will use the simple temporal network from Fig¬ 
ure [^ also for the illustration of some other algorithms 
because it allows the users to manually check the pre¬ 
sented results. 


VI. TEMPORAL CO-OCCURRENCE NETWORK S 

Let the binary matrix A = [oep] describe a two-mode 
network on the set of events E and the set of of partici¬ 
pants P: 


— 


1 p participated in the event e 
0 otherwise 


The function d : E ^ T assigns to each event e the date 
(i(e) when it happened. T = [first, last]. Using these 
data we can construct two temporal affiliation matrices: 


instantaneous Ai = 


^ep 


, where 


^ep 


[{d{e),d{e) + 1,1)] Uep = 1 
\ ] otherwise 


• cumulative Ac = [acep], where 

[(d(e), (asf-I-1,1)] ttep = 1 


^ep 


otherwise 


Using the multiplication of temporal matrices over the 
combinatorial semiring we get the corresponding instan¬ 
taneous and cumulative co-occurrence matrices 


Ci = AU Ai 


and 


Cc = Ac' Ac 


A typical example of such a matrix is the papers author¬ 
ship matrix where E is the set of papers, P is the set of 
authors and d is the publication year ||3]| . 

The triple {s,f,v) in a temporal quantity cipq tells 
that in the time interval [s, /) there were v events in 
which both p and q took part. 

The triple {s,f,v) in a temporal quantity ccpq tells 
that in the time interval [s, f) there were in total v accu¬ 
mulated events in which both p and q took part. 

The diagonal matrix entries cipp and ccpp contain the 
temporal quantities counting the number of events in the 
time intervals in which the participant p took part. 

For example, in a data set on the stem cell research 
during 1997-2012 in Spain collected by Gisela Cantos- 
Mateos 0 we get from the basic two-mode network, 
where E is the set of papers and P is the set of insti¬ 
tutions, for selected two institutions (HCL/B = Univer¬ 
sity Flospital Clinic de Barcelona, Barcelona and IDI/B 
= Institut d’lnvestigacions Biomediques August Pi i 
Sunyer, Barcelona) the collaboration temporal quanti¬ 
ties presented in Table jH] 
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FIG. 4. Intensity of violent activities of police, fascists and all. 


The first column in the table contains the yearly col¬ 
laboration (co-authorship) data and the second column 
contains the cumulative collaboration data. Let’s read 
the table; 

ci[lDI/B, HCL/B](2005, 2006) = 3 — in the year 
2005 researchers from both institutions published 3 joint 
papers; 

ci[lDI/B, HCL/b](2011, 2013) = 18 — in the years 
2011 and 2012 researchers from both institutions pub¬ 
lished 18 joint papers each year; 
ci[HCL/B,HCL/B](2010, 2011) = 78 — in the year 
2010 researchers from the institution HCL/B published 
78 papers; 

cc[IDI/B,HCL/b](2008, 2009) = 16 — till the year 
2008 (included) researchers from both institutions pub¬ 
lished 16 joint papers. 

Note that the violence network from Section|IV]is es¬ 
sentially a co-occurrence network that could be obtained 
from the more primitive instantaneous two-mode net¬ 
work about violent actions reported in journal articles 
and the involved political actors. 


VII. CLUSTERING COEFFICIENTS 

Let us assume that the network Af is based on a sim¬ 
ple directed graph Q = (V, .4) without loops. From 
a simple undirected graph we obtain the corresponding 
simple directed graph by replacing each edge with a pair 
of opposite arcs. In such a graph the clustering coeffi¬ 
cient, C{v), of the node v is defined as the proportion 
between the number of realized arcs among the node’s 
neighbors and the number of all possible arcs among the 
node’s neighbors N{v), that is 

where k is the number of neighbors of the node v. For a 
node V without neighbors or with a single neighbor we 
set C{v) = 0. 

The clustering coefficient measures a local density of 
the node’s neighborhood. A problem with its applica¬ 
tions in network analysis is that the identified densest 
neighborhoods are mostly very small. For this reason 
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FIG. 5. First example network. All unlabeled links have a 
value of [(1, 9,1)]. 


To count the number of realized arcs among the 
node’s neighbors we use the observation that each arc 
forms a triangle with links from its end-nodes to the 
node v\ and that the number of triangles in a simple 
undirected graph can be obtained as the diagonal value 
in the third power of the graph matrix (over the combi¬ 
natorial semiring). 



AAA 

TTT 


AAT 

ATT 


TAT 

ATA 



TTT 

AAA 


ATT 

AAT 


ATA 

TAT 


TABLE II. Temporal collaboration. 


ci [IDI/B,HCL/B] CC [IDI/B,HCL/B] 


1 

(2003, 

2004, 

1) 

(2003, 

2004, 

1) 

2 

(2004, 

2005, 

2) 

(2004, 

2005, 

3) 

3 

(2005, 

2006, 

3) 

(2005, 

2006, 

6) 

4 

(2006, 

2007, 

2) 

(2006, 

2007, 

8) 

5 

(2007, 

2008, 

1) 

(2007, 

2008, 

9) 

6 

(2008, 

2009, 

7) 

(2008, 

2009, 

16) 

7 

(2009, 

2010, 

6) 

(2009, 

2010, 

22) 

8 

(2010, 

2011, 

7) 

(2010, 

2011, 

29) 

9 

(2011, 

2013, 

18) 

(2011, 

2012, 

47) 

10 




(2012, 

2013, 

65) 


ci[HCL/B,HCL/B] 

CC [HCL/B,HCL/B] 

1 

(1997, 

1998, 

2) 

(1997, 

1998, 

2) 

2 

(1998, 

1999, 

5) 

(1998, 

1999, 

7) 

3 

(1999, 

2000, 

8) 

(1999, 

2000, 

15) 

4 

(2000, 

2001, 

7) 

(2000, 

2001, 

22) 

5 

(2001, 

2002, 

5) 

(2001, 

2002, 

27) 

6 

(2002, 

2003, 

6) 

(2002, 

2003, 

33) 

7 

(2003, 

2004, 

14) 

(2003, 

2004, 

47) 

8 

(2004, 

2005, 

20) 

(2004, 

2005, 

67) 

9 

(2005, 

2006, 

10) 

(2005, 

2006, 

77) 

10 

(2006, 

2007, 

14) 

(2006, 

2007, 

91) 

11 

(2007, 

2008, 

20) 

(2007, 

2008, 

111) 

12 

(2008, 

2009, 

28) 

(2008, 

2009, 

139) 

13 

(2009, 

2010, 

56) 

(2009, 

2010, 

195) 

14 

(2010, 

2011, 

78) 

(2010, 

2011, 

273) 

15 

(2011, 

2012, 

84) 

(2011, 

2012, 

357) 

16 

(2012, 

2013, 

112) 

(2012, 

2013, 

469) 


we provided in Pajek the corrected clustering coeffi¬ 
cient, C'{v), 


C'iv) 


l-4(jV(t^))| 

A(fc- 1) 


where A is the maximum number of neighbors in the 
network. 



TAA 

TTA 



TTA 

TAA 


FIG. 6. Counting triangles. 


For simple directed graphs the counting of triangles 
is slightly more complicated. Let us denote T = 
and S = A -f T. From Figure]^ we see that each tri¬ 
angle (determined with a link opposite to the dark node) 
appears exactly once in 

AAA+AAT+TAT+TAA = AAS+TAS = SAS. 

This gives us a simple way to count the triangles which 
is used in Algorithm The function nRows{A) re¬ 
turns the size (number of rows) of matrix A. The 
function VecConstin, v) constructs a vector of size n 
filled with the value v. The function MatBin{A) trans¬ 
forms all values in the triples in the matrix A to 1. 
The function MatSetDiag{A,c) sets all the diagonal 
entries of the matrix A to the value c. The function 
MatSym{A) makes the transformation S = A 0 T. 
Functions VecSum and VecProd implement a component 
wise composition of temporal vectors: VecSum{a, b) = 
[oi (B hi, i = l,...,n] and VecProd{a,b) = 

[tti Q bi, i = l,...,n]. Similarly Veclnv{a) = 
[invert{ai), i = 1,... ,n] in the combinatorial semir¬ 
ing; where invert{a) = [(s,/, 1/u) for (s,/, u) G 

a]. The function MatProd{A, B) determines the prod¬ 
uct A © B. Since we need only the diagonal val¬ 
ues of the matrix SAS we applied a special function 
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MatProdDiag{A, B) that determines only the diago¬ 
nal vector of the product A © B. Afterward, to get 
the clustering coefficient, we have to normalize the ob¬ 
tained counts. The number of neighbors of the node v 
is determined as its degree in the corresponding undi¬ 
rected temporal skeleton graph (in which an edge e = 
{v : u) exists iff there is at least one arc between the 
nodes v and u). The maximum number of neighbors 
A can be considered either for a selected time point 
{type = 2) or for the complete time window {type = 3). 
Note that to determine the temporal A we used sum¬ 
ming of temporal degrees over the maxmin semiring 


(M, max, mm, —oo, oo). 

The time complexity of Algorithm[^is 0(n^ • L). 

In Table III and Table the ordinary and the cor¬ 
rected clustering coefficients are presented for the ex¬ 
ample network from Figure and its undirected skele¬ 
ton. 


TABLE III. Clustering coefficients for the first example net¬ 
work. 

1 : [] 

2 : [] 

3 : [] 


4 

[ 

(1, 

3, 

0 

5), (3, 

9, 

0 . 

1667) 

5 

[ 

(1, 

5, 

0 

1667) , 

(5, 

9, 

0 

5) 

6 

[ 

(1, 

9, 

0 

5) ] 





7 

[ 

(1, 

5, 

0 

25), (5 

, 9 

0 

.5) 

] 

8 

[ 

(1, 

7, 

0 

4167), 

(7, 

9, 

0 

5) 

9 

[ 

(1, 

7, 

0 

4167), 

(7, 

9, 

0 

5) 

10 

[ 

(1, 

7, 

0 

4167), 

(7, 

9, 

0 

5) 

11 

[ 

(1, 

9, 

0 

5) ] 





12 

[ 

] 








13 

[ 

(2, 

8, 

1 

0) ] 





14 

[ 

(2, 

8, 

1 

0) ] 





15 

[ 

(2, 

8, 

1 

0) ] 






TABLE IV. Corrected clustering coefficients for the skeleton 
of the first example network. 


ALGORITHM 3. Clustering coefficients. 

1: function clusCoef {A, type — 1) 

2: # type — 1 - standard CC 
3: # type = 2 - corrected CC / temporal degMax 
4: # type = 3 - corrected CC / overall degMax 
5: SetSemiring{combinatorial) 

6: n t— nRows{A); ve t— VecConst{n, [(0, oo, —1)]) 

7: B <— MatSetDiag{MatBin{A),0) 

8: St— MatBin{MatSym{B)) 

9: deg t— MatVecMulR{S, VecConst{n, 1)) 

10: if type = 1 then 

11: fac-ir- VecProd{deg,VecSum{deg,ve)) 

12: else 

13: SetSemiring{maxmin); <5 t— 0 

14: for d G deg do <5 t— sum{S, d) 

15: if/ype = 3 then 

16: A t— max([u for (s, /, v) G 5]) 

17: (5 t- [(0,oo, A)] 

18: SetSemiring{combinatorial) 

19: degm t— VecSum{deg,ve)',fac t— 0 

20: for d G degm do fac.append(pwd{5, d)) 

21: tri t— MatProdDiag{MatProd{S, B), S) 

22: return VecProd{VecInv{fac),tri) 


1 : [] 
2 : [] 
3 : [] 


4 

[ 

(1, 

3, 

0 

5), (3, 

9, 

0 . 

25) 

] 

5 

[ 

(1, 

5, 

0 

25), (5 

, 9, 

0 

.5) 

] 

6 

[ 

(1, 

9, 

0 

5) ] 





7 

[ 

(1, 

5, 

0 

5), (5, 

7, 

0 . 

75) 

f 



(7, 

9, 

1 

0) ] 





8 

[ 

(1, 

7, 

0 

8333), 

(7, 

9, 

1. 

0) 

9 

[ 

(1, 

7, 

0 

8333), 

(7, 

9, 

1. 

0) 

10 

[ 

(1, 

7, 

0 

8333), 

(7, 

9, 

1. 

0) 

11 

[ 

(1, 

7, 

0 

75), (7 

, 9, 

1 

.0) 

] 

12 

[ 

] 








13 

[ 

(2, 

8, 

0 

5) ] 





14 

[ 

(2, 

8, 

0 

5) ] 





15 

[ 

(2, 

8, 

0 

5) ] 






In some semirings different closures can exist. For com¬ 
puting the matrix closure we can apply the Fletcher’s al¬ 
gorithm Qo). The entry Cuv in the matrix C = A* is 
equal to the sum of values of all walks from the node u 
to the node v. In most of the semirings, except the com¬ 
binatorial, for which we are interested in determining 
the closures, also the absorption law holds 


1 © a = 1, for all a G A. 


VIII. CLOSURES IN TEMPORAL NETWORK S 

When the basic semiring (A, ©, ©, 0,1) is closed - 
an unary closure operation * with the property 

a* = 1 © a © a* = 1 © a* 0 a, for all a C A 

is defined in it - this property can be extended also to 
the corresponding matrix semiring. When it exists, a 
standard closure is obtained as 

OO 

a* = 0a\ 

i=0 


In these semirings a* = 1, for all a G A, and there¬ 
fore the Fletcher’s algorithm can be simplified and per¬ 
formed in place as implemented in Algorithm]^ 

For a temporal quantity a over a closed semiring it 
holds Ta* = T. 

The time complexity of Algorithm|^is 0{n^ ■ L). 

IX. TEMPORAL NODE PARTITIONS 

In the previous sections, the nodes of temporal net¬ 
works were considered as being present all the time. 
We can describe the presence of nodes through time 
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ALGORITHM 4. Closure of a temporal matrix over an ab¬ 
sorptive semiring. 

1: function MatClosure{R, strict = False) 

2: n nRows{R) 

3: C ^R 

4: for fc G 1 : n do 

5: for n G 1 : n do 

6: for G 1 : n do 

7: C[m, u] •<— 

8: ™m(C[u, u],proc/(C'[M, fc], C[k, w])) 

9: if -^strict then C[k, k] sum{l, C[k, k]) 

10: return C 

using a temporal binary (single valued) node partition 
T : V ^ A^iT), 

T{u) = forwGV 

specifying that a node u is present in time intervals 

[si,fi),i = 

The node partition Tmiu determined from the tempo¬ 
ral network links by 

u binary (ai), 

l£C:uGext{l) 

for u G V, is the smallest temporal partition of nodes 
that satishes the consistency condition from Section [n| 
The term ext(Z) denotes the set of endnodes of the link 
I, ai is the temporal quantity assigned to the link I, and 
the function binary sets all values in a given temporal 
quantity to 1. In the library TQ the partition TMin can 
be computed using the function minTime. 

A temporal node partition q can also be used to ex¬ 
tract a corresponding subnetwork from the given tem¬ 
poral network described with a matrix A. The subnet¬ 
work contains only the nodes active in the partition q 
and the active links satisfying the consistency condition 
with respect to q. 

To formalize the described procedure we first define 
the procedure extract{p, a) = b, where p is a binary tem¬ 
poral quantity and a is a temporal quantity, as 

t G TpHTa 

'' ' [ i€ otherwise 

Let B be a temporal matrix describing the links of the 
subnetwork determined by the partition q. Its entries for 
l{u, v) G C are determined by 

bi = extract{q(u) H q{v),ai). 

In TQ this operation is implemented as a procedure 
MatExtract{c{, A). 


TABLE V. Temporal input reachability degrees for the first 
example network. 


1 

[ (1, 

9, 

3) ] 






2 

[ (1, 

9, 

3) ] 






3 

[] 








4 

[ (1, 

3, 

3) , 

(3, 9, 

6) ] 




5 

[ (1, 

3, 

3) , 

(3, 9, 

6) ] 




6 

[ (1, 

3, 

3) , 

(3, 9, 

6) ] 




7 

[ (1, 

3, 

3) , 

(3, 5, 

6) , 

(7, 

9, 

5) ] 

8 

[ (1, 

3, 

8) , 

(3, 5, 

11) 

(5, 

9, 

5) 

9 

[ (1, 

3, 

8) , 

(3, 5, 

11) 

(5, 

9, 

5) 

10 

[ (1, 

3, 

8) , 

(3, 5, 

11) 

(5, 

9, 

5) 

11 

[ (1, 

3, 

8) , 

(3, 5, 

11) 

(5, 

9, 

5) 

12 

[] 








13 

[ (2, 

8, 

3) ] 






14 

[ (2, 

8, 

3) ] 






15 

[ (2, 

8, 

3) ] 







X. TEMPORAL REACHABILITY AND WEAK AND 
STRONG CONNECTIVITY 

For a temporal network represented with the cor¬ 
responding binary matrix A its transitive closure A* 
(over the reachability semirings based on the semiring 
({0,1}, V, A, 0,1)) determines its reachability relation 
matrix. We obtain its weak connectivity temporal ma¬ 
trix W as 

W= (AU A^)* 

and its strong connectivity temporal matrix S as 
S = A* n (A*)^. 

The use of the strict transitive closure instead of a tran¬ 
sitive closure in these relations preserves the inactivity 
value 0 on the diagonal for all isolated nodes. 

A. Reachability degrees 

Let R = A = A © A* be the strict reachability re¬ 
lation of a given network. Then the temporal vectors 
inReach = /nDeg(R) and outReach — outDeg(R) con¬ 
tain temporal quantities counting the number of nodes; 
from which a given node v is reachable {inReach[v]) 
/ which are reachable from the node v (outReach[v]). 
The results for our example network are presented in 
Table |v] For example, 8 nodes {4, 5,6, 7,8, 9,10,11} 
are reachable from node 6 in the time interval [1, 5), and 
3 nodes {4, 5,6} are reachable in the time interval [5, 9). 

B. Temporal weak connectivity 

The function weakConnMat{A) for a given temporal 
network matrix A determines the corresponding tempo¬ 
ral weak connectivity matrix W. Every time slice 
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t G T, of the matrix W is an equivalence relation that 
can be compactly described with the cori'esponding par¬ 
tition. 

To transform the temporal equivalence matrix E into 
the corresponding temporal partition p we use the fact 
that on a given time interval equivalent (in our case 
weakly connected) nodes get the same value on this in¬ 
terval in the product of the matrix E with a vector com¬ 
puted over the combinatorial semiring (N,-f, •, 0,1). 
We take for the vector values randomly shuffled inte¬ 
gers from the interval 1 : n. With a very high prob¬ 
ability the values belonging to different equivalence 
classes are different. This is implemented as a procedure 
eqMat2Part{E) (see Algorithm|^. Maybe in the future 
implementations we shall add a loop with the check of 
the injectivity of this mapping. The classes of the ob¬ 
tained temporal partition are Anally renumbered with 
consecutive numbers using the function renumPart{p) 
(see Algorithm]^. The variable C in the description 
of the function renumPart is a dictionary (data struc¬ 
ture). 

For our first example network we obtain the tempo¬ 
ral weak partition presented on the left hand side of Ta¬ 
ble IVj 

ALGORITHM 5. Transform temporal equivalence relation 
into partition. 

I: function eqMatlPartiE) 

2: SetSemiring{combinatorial) 

3: u shuffle{\\(f), oo, i + 1)] for i G 1 : nRows{E)]) 

4: p MalVecMulR{E, u) 

5: return renumPart(p) 


ALGORITHM 6. Renumber the classes of a partition. 

1: function renumPart(p) 

2: C^{};g=[] 

3: for a G p do 

4: r.(-[] 

5: for (Sa,/a, Ca) € a do 

6: if Ca ^ C then C[ca] t— 1 -f length{C) 

7: r.append{{sa, fa,C[Ca])) 

8: q.append{r) 

9: return q 


TABLE VI. Temporal weak and strong connectivity partitions 
for the first example network. 


Weak partition 


1 

[ (1, 

3, 

1) , 

(3, 5, 

2) , 

(5, 

9, 

3) 

2 

[ (1, 

3, 

1) , 

(3, 5, 

2) , 

(5, 

9, 

3) 

3 

[ (1, 

3, 

1) , 

(3, 5, 

2) , 

(5, 

9, 

3) 

4 

[ (1, 

3, 

4) , 

(3, 5, 

2) , 

(5, 

9, 

3) 

5 

[ (1, 

3, 

4) , 

(3, 5, 

2) , 

(5, 

9, 

3) 

6 

[ (1, 

3, 

4) , 

(3, 5, 

2) , 

(5, 

9, 

3) 

7 

[ (1, 

3, 

4) , 

(3, 5, 

2) , 

(5, 

9, 

5) 

8 

[ (1, 

3, 

4) , 

(3, 5, 

2) , 

(5, 

9, 

5) 

9 

[ (1, 

3, 

4) , 

(3, 5, 

2) , 

(5, 

9, 

5) 

10 

[ (1, 

3, 

4) , 

(3, 5, 

2) , 

(5, 

9, 

5) 

11 

[ (1, 

3, 

4) , 

(3, 5, 

2) , 

(5, 

9, 

5) 

12 

[] 








13 

[ (2, 

8, 

6) ] 






14 

[ (2, 

8, 

6) ] 






15 

[ (2, 

8, 

6) ] 






Strong partition 






1 

[ (1, 

9, 

1) ] 






2 

[ (1, 

9, 

1) ] 






3 

[] 








4 

[ (1, 

9, 

2) ] 






5 

[ (1, 

9, 

2) ] 






6 

[ (1, 

9, 

2) ] 






7 

[ (7, 

9, 

3) ] 






8 

[ (1, 

7, 

4) , 

(7, 9, 

3) ] 




9 

[ (1, 

7, 

4) , 

(7, 9, 

3) ] 




10 

[ (1, 

7, 

4) , 

(7, 9, 

3) ] 




11 

[ (1, 

7, 

4) , 

(7, 9, 

3) ] 




12 

[] 








13 

[ (2, 

8, 

5) ] 






14 

[ (2, 

8, 

5) ] 






15 

[ (2, 

8, 

5) ] 







The time complexity of algorithms for temporal weak 
and strong connectivity partitions is 0{'n? ■ L). 

For our first example network we obtain the tempo¬ 
ral strong partition presented on the right hand side of 
Table VI In the library TQ both matrices and partitions 
are based on the strict transitive closure. 


XI. TEMPORAL CLOSENESS AND BETWEENNESS 

Closeness and betweenness are among the traditional 
social network analysis indices measuring the impor¬ 
tance of nodes ca. They are somehow problematic 
when applied to non (strongly) connected graphs. In this 
section we will not consider these questions. We will 
only show how to compute them for non-problematic 
temporal graphs. 


C. Temporal strong connectivity 


The procedure strongConnMat{A) for a given tem¬ 
poral network matrix A determines the corresponding 
temporal strong connectivity matrix S. To determine the 
intersection of temporal network binary matrices A and 
B we use the function MatInter{A, B). Again, to get 
the strong connectivity partition we have to apply the 
function eqMatIPart to the strong connectivity matrix. 


A. Temporal closeness 


The output closeness of the node v is defined as 

It \ 

ocl[v) = 


^ ^ dvu 
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FIG. 7. Second example network. All unlabeled arcs have the 
value [(1,9,1)]. 


To determine the closeness we first need to compute 
the matrix D = [duv] of geodetic distances duv between 
the nodes u and v. It can be obtained as a closure of 
the network matrix A over the shortest paths semiring 
(MJ, min, +, 00 , 0). Note that the values in the matrix 
A can be any nonnegative real numbers. 

In Figure]^ we present our second example temporal 
network which is an extended version of the example 
given in Figure 3 from IT] . 

Because a complete strict closure matrix D is too 
large to be listed we present only some of its selected 
entries: 


D[3,1] 

= [ (3, 

7, 

3) , 

(7, 9, 

5) 

D[4,6] 

= [ (1, 

4, 

1) , 

(4, 6, 

5) 

D[6,3] 

= [ (3, 

5, 

6) , 

(5, 9, 

4) 

D[7,6] 

= [ (1, 

9, 

4) ] 




To compute the vector of closeness coefficients of 
nodes we have to sum the temporal distances to other 
nodes over the combinatorial semiring. See Algo¬ 
rithm]^ While summing we replace gaps (inactivity in¬ 
tervals inside T) with time intervals with the value infin¬ 
ity, using the procedure The time complexity 

of Algorithm]^ is 0(n^ ■ L). 

The temporal closeness coefficients for our second 


example network are given in Table VII 


ALGORITHM 7. Temporal closeness. 

1: function ctojenewf A, type = 2) 

2: # type'. 1 - output, 2 - all, 3 - input 
3: St— startTime{A)', f ^ finishTime{A) 

4: n t— nRows{A) 

5: SetSemiring{path) 

6: D <— MatClosure{A, strict — True) 

1: SetSemiring{combinatorial) 

8: fc t— (2 — \type — 2|) • (n — l);/flc t— [(0, oo, k)] 

9: for u € 1 : n do 

10: d t- 0 

11: for u G 1 : n do 

12: itu ^ V then 

13: if type < 3 then 

14: d sum(d,fillGaps{D[v,u\, s, f)) 

15: if type > 1 then 

16: d<^sum(d,fillGaps{D[u,v],s,f)) 

17: cl[v\ t— prod(fac,mvert{d)) 

18: return cl 


TABLE VII. Output closeness for the second example net¬ 
work. 


1 ; 

[ (1, 

9, 

0.4375)] 




2 : 

[(1, 3, 
(5, 9, 

0.0000) , 
0.5833)] 

(3, 

5, 

0.4375) 

3 ; 

[(1, 3, 
(7, 9, 

0.0000), 
0.3889) ] 

(3, 

7, 

0.4375) 

4 ; 

[ (1, 

3, 

0.0000), 

(3, 

4, 

0.4375) 


(4, 6, 
(7, 9, 

0.3500), 
0.3500) ] 

(6, 

7, 

0.4375) 

5 : 

[(1, 3, 
(7, 9, 

0.0000), 
0.3500)] 

(3, 

7, 

0.4375) 

6 : 

[(1, 3, 
(5, 9, 

0.0000) , 
0.3500)] 

(3, 

5, 

0.2917) 

7 : 

[(1, 3, 
(7, 9, 

0.0000) , 
0.3500)] 

(3, 

7, 

0.4375) 

8 ; 

[(1, 3, 
(5, 9, 

0.0000) , 
0.4375)] 

(3, 

5, 

0.3500) 


where du,v is the length of u-v geodesics. Then it is also 
easy to determine the quantity nu,w{v): 


nu,w{v) 


ttu^V ‘ dy.^y “t“ dy - du.W 

0 otherwise 


B. Temporal betweenness 

The betweenness of a node v is defined as 

. / N ^ nu,w{v) 

(n-l)(n-2) nu,w 

|{„“C™}I=3 

where nu,w is the number of u-w geodesics (shortest 
paths) and is the number of u-w geodesics pass¬ 

ing through the node v. 

Suppose that we know the matrix 

C {{dy^y ^ Tly^y)^ 


This gives the following scheme of procedure for com¬ 
puting the nontemporal betweenness coefficients b 

1: compute C 
2 : for u e V do 

3: r ^ 0 

4: tor u €V,w gV do 

5: if n[u,w] ^ 0 A |{t;,u,w}| = 3 A 

d[u, w] = d[u, u] + d[v, w] then 
6: rr + n[u,v] ■ n[v,w]/n[u,w] 

7: b[v] G- rj ((n — 1) • (n — 2)) 

In HI it is shown that the matrix C can be obtained 
by computing the closure of the network matrix over 












13 


_2 

the geodetic semiring (N , 0 , ©, (oo, 0), (0,1)), where 
N = N U {c»} and we dehne addition © with 


i i a < b 

i+j a = b) 

j a>b 


and multiplication © with; 

{a,i)Q{b,j) = {a + b,i-j). 


To compute the geodetic closure we hrst transform 
the network temporal adjacency matrix A to a matrix 
G = [(d, ri)u,v] which has for entries pairs dehned by 


(d, Tl)u^y (f) 


(1,1) 31 G £ : l{u,v) A t G T{1) 

otherwise 


where d is the length of a geodesic and n is the num¬ 
ber of geodesics from uto v. In temporal networks the 
distance d and the counter n are temporal quantities. 

The presented scheme adapted for computing the 
temporal betweenness vector is implemented in TQ as 
the function betweenness{A). First we compute its strict 
geodetic closure C over the geodetic semiring. We 
present only some of its selected entries for our second 
example network; 


C[l,7] 

= 1(1, 9, 

(3, 4))] 



C[2,2] 

= [ (1, 3, 

(4, 4) ) , 

(3, 4, 

(4, 6)), 


(4, 5, 

(4, 5)), 

(5, 9, 

(2, 1))] 

C[4,6] 

= [ (1, 4, 

(1, 1) ) , 

(4, 6, 

(5, 3)), 


(6, 9, 

(1, 1) ) ] 



C[5,5] 

= 1(1, 9, 

(1, 1) ) ] 



C[6,3] 

= 1(3, 5, 

(6, 2) ) , 

(5, 9, 

(4, 1) ) ] 

C[7,6] 

= 1(1, 3, 

(4, 2) ) , 

(3, 4, 

(4, 6)), 


(4, 6, 

(4, 3)), 

(6, 7, 

(4, 6)), 


(7, 9, 

(4, 2) ) ] 




For example, the value C[4,6] reflects the facts that 
an arc exists from node 4 to node 6 in time intervals 
[1,4) and [6,9); and in the time interval [4,6) they are 
connected with 3 geodesics of length 5; (4, 7,8, 2, 5,6), 
(4,7,1,3,5,6), (4,7,1,2,5,6). 

We continue and using the combinatorial semiring we 
compute the temporal betweenness vector b. The speci¬ 
ficity of temporal quantities d[u, z;] and n[u, z;] is consid¬ 
ered in the auxiliary function between that implements 
the temporal version of the statement 

if d[u, w] = d[u, w] + d[v, w] then 
r ^ r + n[u, u] • n[v, w\/n[u, w] 


from the basic betweenness algorithm. Again we apply 
the merging scheme. The time complexity of the proce¬ 
dure betweenness is 0{n^ ■ L). 


The temporal betweenness coefficients for our second 


example network are presented in Table VIII 


TABLE VIII. Betweenness for the second example network. 

1 : 

[ (3, 

4, 

0. 

,2500) , 

(4, 

6, 

0 , 

.2754), 


(6, 

7, 

0. 

.2500), 

(7, 

9, 

0 , 

.1429)] 

2 : 

[ (1, 

3, 

0. 

.3452), 

(3, 

4, 

0 , 

.4048), 


(4, 

6, 

0. 

.4187) , 

(6, 

7, 

0 , 

. 4048), 


(7, 

9, 

0. 

. 6071) ] 





3 : 

[ (1, 

3, 

0. 

,0595) , 

(3, 

4, 

0 , 

.0952), 


(4, 

6, 

0. 

.1052) , 

(6, 

7, 

0 , 

.0952), 


(7, 

9, 

0. 

,0595) ] 





4 : 

[ (1, 

3, 

0. 

.1667), 

(3, 

4, 

0 , 

.2500), 


(4, 

5, 

0. 

.1762), 

(5, 

6, 

0 , 

. 1048), 


(6, 

9, 

0. 

,1786)] 





5 : 

[ (1, 

3, 

0. 

.1667), 

(3, 

4, 

0 , 

.2500), 


(4, 

5, 

0. 

,3476), 

(5, 

6, 

0 , 

.2762), 


(6, 

9, 

0. 

.1786)] 





6 : 

[ (1, 

3, 

0. 

,1190), 

(3, 

4, 

0 , 

.0952), 


(4, 

6, 

0. 

.0544) , 

(6, 

7, 

0 , 

.0952), 


(7, 

9, 

0. 

,1786)] 





7 : 

[ (1, 

3, 

0. 

.1190), 

(3, 

4, 

0 , 

. 4048), 


(4, 

5, 

0. 

.4694), 

(5, 

6, 

0 , 

.3266), 


(6, 

7, 

0. 

,2619), 

(7, 

9, 

0 , 

. 1786)] 

8 : 

[ (1, 

3, 

0. 

.3095), 

(3, 

4, 

0 , 

.2500), 


(4, 

6, 

0. 

,2484) , 

(6, 

7, 

0 , 

.2500), 


(7, 

9, 

0. 

.5238) ] 






XII. TEMPORAL PATHFINDER 

The Pathhnder algorithm was proposed in the eight¬ 
ies EtI l28l for the simplihcation of weighted networks 
- it removes from the network all links that do not 
satisfy the (generalized) triangle inequality - if for a 
weighted link there exists a shorter path connecting its 
endnodes then the link is removed. The basic idea of the 
Pathhnder algorithm is simple. It produces a network 
PFnet(W,r,(z) = {V,Cpf) determined by the follow¬ 
ing scheme of procedure 

1: compute 
2 : CpF 'G- 0 ; 

3: for e{u, u) G £ do 

4: if [y^ = W[m, u] then 

5: £pF G- £pF U {e} 


where W is a network dissimilarity matrix and 'W^i') = 
0^^^ W* = (1 © W)"? is the matrix of the values of all 
walks of length at most q computed over the Pathfinder 
semiring ©, 0 , oo, 0 ) with a\nb = s/ajAfifir 
a © 6 = min(a, b). The value of Wuv{q) in the matrix 
'W^i) is equal to the value of all walks of length at most 
q from the node u to the node v. 

The scheme of Pathhnder is implemented as the func¬ 
tion pathFinder. The temporal version of the statement 

if W(9)[it,z;] = W[u,v] then £pf := £pf U {e} 
is implemented in the function PFcheck using the merg¬ 
ing scheme. 

The function MatPower{A, k) computes the fc-th 
power of the matrix A. 

The time complexity of Algorithm [sj+j^ is 0{L ■ rA ■ 
logq) US). 
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ALGORITHM 8. Temporal PathFinder. 

I: function pathFinder(W, r — 1, q = oo) 

2: n ■<— nRows{W)\ SetSemiringipathfinder, r, q) 

3: if g > n then Z MatClosure{W) 

4: else Z -ir- MatPower{MatSetDiag{W, 1), q) 

5: for u ^ 1 : n,v ^ 1 : n do 

6: PF[u, n] •<— PFcheck{W[u, n], Z[u, ?;]) 

7: return PF 


r = 1 a link e is removed if there exists a path, con¬ 
necting its initial node to its terminal node, with the 
value (sum of link values) smaller than the value of 
the link e. The arc (1,2) is removed because 3 = 
u(l, 2) > u(l, 3) -I- u(3, 2) = 2. The arc (1, 6) is re¬ 
moved in the time interval [5,9) because in this interval 
5 = u(l,6) > u(l,3)-ft;(3,4)-Fu(4,5)-fu(5,6) =4. 


ALGORITHM 9. Temporal PathFinder merge operation. 
1: function PFcheck(a,b) 

2: if length{a) — 0 then return a 

3: if length{b) — 0 then return a 

4: C t— [ ] 

5: {Sa,fa,Va) ^ get{a)-, {sb,fb,Vb) <-get{b) 

6: while {sa < oo) V {sb < oo) do 

7: if fa < Sb then (sa, fa,Va) get{a) 

8: else if fb < Sa then (sb, fb, Vb) t— get{b) 

9: else 

10: Sc t- max(sa, Sb); fc fb) 

11: if Vb = Va then c.append{{sc, fc, Va)) 

12: if fc = fa then (s^, fa, Va) ^ get{a) 

13: if fc = fb then (sb, fb,Vb)get{b) 

14: return standardise) 



FIG. 8. Pathfinder example. 


The bottom network in Figure presents the 
Pathfinder skeleton PFnet(A/^, 1, oo) of a network N 
presented in the top part of the same figure. Because 


XTTT. SEPTEMBER 1ITH REUTERS TERROR NEWS 


The Reuters terror news network was obtained from 
the CRA (Centering Resonance Analysis) networks pro¬ 
duced by Steve Corman and Kevin Dooley at Arizona 
State University. The network is based on all the stories 
released during 66 consecutive days by the news agency 
Reuters concerning the September 11 attack on the U.S., 
beginning at 9:00 AM EST 9/11/01. The nodes of this 
network are important words (terms). There is an edge 
between two words iff they appear in the same utterance 
(for details see the paper jO)). The weight of an edge 
is its frequency. The network has n = 13332 nodes 
(different words in the news) and m = 243447 edges, 
50859 with value larger than 1. There are no loops in 
the network. 

The Reuters terror news network was used as a 
case network for the Viszards visualization session on 
the Sunbelt XXII International Sunbelt Social Network 
Conference, New Orleans, USA, 13-17. February 2002. 

We transformed the Pajek version of the network into 
the lanus format used in TQ. To identify important terms 
we computed their aggregated frequencies and extracted 
the subnetwork of the 50 most frequently used (during 


66 days) nodes. They are listed in Table IX 


Trying to draw this subnetwork it turns out to be al¬ 
most a complete graph. To obtain something readable 
we removed all temporal edges with a value smaller than 
10. The corresponding underlying graph is presented in 
Figure]^ The isolated nodes were removed. 

For each of the 50 nodes we determined its temporal 
activity and drew it. By visual inspection we identified 6 
typical activity patterns - types of terms (see Figure [TOji. 
For all charts in the figure the displayed values are in the 
interval [0, 200] - the largest activity value for the term 
Wednesday is larger than 200. 

The primary terms are the terms with a very high fre¬ 
quency of appearance in the first week after September 
11th and smaller, slowly declining values in the follow¬ 
ing period. The representative of this group in Figure [T0| 
is hijack and other members are: airport, american, at¬ 
tack, city, day, flight, nation. New York, official. Pen¬ 
tagon, people, plane, police, president Bush, security, 
tower. United States, Washington, world. World Trade 
center. These are the terms describing the event. 

The secondary terms are a reaction to the event. 
There are no big changes in their values. We identified 
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FIG. 9. September 11th. 


three subgroups: a) slowly declining represented with 
bin Laden (country, foreign, government, military, min¬ 
ister, new, Pakistan, tell, terrorism, terrorist, time, war, 
week); b) stationary represented with taliban (afghan, 
Afghanistan, force, group, leader); and c) occasional 
with several peaks, represented with bomb (air, build¬ 
ing, office, strike, worker). 

There are three special patterns - two periodic 
Wednesday and Tuesday; and one episodic anthrax. 

To consider in a measure of importance of the node 
u G V also the node’s position in the network we con¬ 
structed the attraction coefficient att(M). 

Let A = [quv] be a network matrix of temporal quan¬ 
tities with positive real values. We define the node ac¬ 
tivity act(u) as (see Section IV i 


act(M) = act({ri}, V \ {u}) = ^ a^. 

t)GV\{u} 

Then the attraction of the node u is dehned as 


att(w) = X H 


A act(u) 


Note that the fraction is measuring the proportion 
of the activity of the node v that is shared with the node 
u. 


From 0 < < 1 and deg(u) = 0 ^ = 0 it 

follows that 


E 


^vu 

act(u) 


< deg(u) < A 


where A denotes the maximum degree. Therefore we 
have 0 < att(it) < 1, for all u €V. 

The maximum possible attraction value 1 is attained 
exactly for nodes: a) in an undirected network: that are 
the root of a star; b) in a directed network: that are the 
only out-neighbors of their in-neighbors - the root of a 
directed in-star. 

We computed the temporal attraction and the corre¬ 
sponding aggregated attraction values for all the nodes 
in our network. We selected 30 nodes with the largest 
aggregated attraction values. They are listed in Table 
Again we visually explored them. In Figure [m we 
present temporal attraction coefficients for the 6 selected 
terms. For all charts in the figure the displayed attraction 
values are in the interval [0,0.2]. 

Comparing on the common terms (taliban, bomb, an¬ 
thrax) the activity charts in Figure [TO] w ith the cor¬ 
responding attraction charts in Figure |1 f] we see that 
they are “correlated” (obviously act(a; t) = 0 implies 
att(a; t) = 0), but different in details. 

For example, the terms taliban and bomb have small 
attraction values at the beginning of the time window - 
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hijack : 




I ti haJaJ] n .- ^ — 


bin Laden : 



Wednesday : 



FIG. 10. Types of activity. 


TABLE IX. 50 most frequent terms in the Terror news net¬ 
work. 


n 

term 

Efreq 

n 

term 

Efreq 

1 

united_states 

15000 

26 

terrorism 

2212 

2 

attack 

10348 

27 

day 

2128 

3 

taliban 

6266 

28 

week 

2017 

4 

people 

5286 

29 

worker 

1983 

5 

afghanistan 

5176 

30 

office 

1967 

6 

binjaden 

4885 

31 

group 

1966 

7 

new_york 

4832 

32 

air 

1962 

8 

pres_bush 

4506 

33 

minister 

1919 

9 

Washington 

4047 

34 

time 

1898 

10 

official 

3902 

35 

hijack 

1884 

11 

anthrax 

3563 

36 

strike 

1818 

12 

military 

3394 

37 

afghan 

1775 

13 

plane 

3078 

38 

flight 

1775 

14 

world_trade_ctr 

3006 

39 

tell 

1746 

15 

security 

2906 

40 

terrorist 

1745 

16 

american 

2825 

41 

airport 

1741 

17 

country 

2794 

42 

Pakistan 

1714 

18 

city 

2689 

43 

tower 

1685 

19 

war 

2679 

44 

bomb 

1674 

20 

tuesday 

2635 

45 

new 

1650 

21 

pentagon 

2620 

46 

buildng 

1634 

22 

force 

2516 

47 

Wednesday 

1593 

23 

government 

2380 

48 

nation 

1589 

24 

leader 

2375 

49 

police 

1587 

25 

world 

2213 

50 

foreign 

1558 


the terms were disguised by the primary terms. On the 
other hand, the terms taliban and Kabul get increased 
attraction towards the end of the time window. 


XIV. CONCLUSIONS 

In the paper we proposed an algebraic approach to 
the “deterministic” analysis of temporal networks based 
on temporal quantities and presented algorithms for the 
temporal variants of basic network analysis measures 
and concepts. We expect that the support for many 
temporal variants of other network analysis notions can 
be developed in similar ways. Our results on tempo¬ 
ral variants of eigen values/vectors based indices (Katz, 
Bonacich, hubs and authorities, page rank) are presented 
in a separate paper ||25]| . 

The proposed approach is an alternative to the tradi¬ 
tional cross sectional approach based on time slices. Its 
main advantages are: 


• the data and the results are expressed using tern- 
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pres Bush : 


TABLE X. 30 most attractive terms in the Terror news net¬ 
work. 





L 


taliban : 





FIG. 11. Attraction patterns. 


n 

term 

Eatt 

n 

term 

Eatt 

1 

united_s tales 

12.216 

16 

war 

2.758 

2 

taliban 

7.096 

17 

force 

2.596 

3 

attack 

7.070 

18 

new_york 

2.590 

4 

afghanistan 

5.142 

19 

government 

2.496 

5 

people 

5.023 

20 

day 

2.338 

6 

binJaden 

4.660 

21 

leader 

2.305 

7 

anthrax 

4.601 

22 

terrorism 

2.202 

8 

pres_busb 

4.374 

23 

time 

2.182 

9 

country 

3.317 

24 

group 

2.072 

10 

Washington 

3.067 

25 

afghan 

2.040 

11 

security 

2.939 

26 

world 

1.995 

12 

american 

2.922 

27 

week 

1.961 

13 

official 

2.831 

28 

Pakistan 

1.943 

14 

city 

2.798 

29 

letter 

1.866 

15 

military 

2.793 

30 

new 

1.851 


poral quantities that are natural descriptions of 
properties changing through time; 

• the user does not need to be careful about the in¬ 
tervals on which the time slices are determined 
- exactly the right intervals are selected by the 
merging (sub)operations. This also improves, 
on average, the efficiency of the proposed algo¬ 
rithms. 

All the described algorithms (and some others) are 
implemented in a Python library TQ (Temporal Quan¬ 
tities) available at the page Q. We started to develop 
a program lanus that will provide a user-friendly (Pajek 
like) access to the capabilities of the TQ library. 

The main goal of the paper was to show; it can be 
done. Therefore we based the current version of the li¬ 
brary TQ on a matrix representation of temporal net¬ 
works as it is presented in the paper. For this repre¬ 
sentation most of the network algorithms have the time 
complexity of 0{n^ ■ L) and the space complexity of 
0{n^ ■ L). This implies that their application is lim¬ 
ited to networks of moderate size (up to some thousands 
of nodes). Large networks are usually sparse. On this 
assumption more efficient algorithms can be developed 
based on a graph (sparse matrix) representation - one of 
the directions of our future research. 

In a description of a temporal network A/^ we can con¬ 
sider also a transition time or latency t G W: t(/, t) is 
equal to the time needed to traverse the link I starting 
at the instant t. Problems considering latency are typi¬ 
cal for operations research but could be important, when 
such data are available, also in social network analysis 
iiiniEiiEaiiii. The analysis of temporal networks 
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considering also the latency seems a much harder task 
- for example, in such temporal networks the strongly 
connected components problem is NP-complete Q • 


The results obtained from temporal procedures are 
relatively large. To identify interesting elements we 
used in the paper the aggregated values and the visual¬ 
ization of selected elements. Additional tools for brows¬ 


ing and presenting the results should be developed. 
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